Modifying molecule-surface scattering by ultrashort laser pulses 
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In recent years it became possible to align molecules in free space using ultrashort laser pulses. 
Here we explore two schemes for controlling molecule-surface scattering process, which are based 
on the laser-induced molecular alignment. In the first scheme, a single ultrashort nonresonant 
laser pulse is applied to a molecular beam hitting the surface. This pulse modifies the angular 
distribution of the incident molecules, and causes the scattered molecules to rotate with a preferred 
sense of rotation (clockwise or counter-clockwise). In the second scheme, two properly delayed laser 
pulses are applied to a molecular beam composed of two chemically close molecular species (isotopes, 
or nuclear spin isomers). As the result of the double pulse excitation, these species are selectively 
scattered to different angles after the collision with the surface. These effects may provide new 
means for the analysis and separation of molecular mixtures. 

PACS numbers: 37.10.Vz,34.35.+a,34.50.Rk,68.49.Df 
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I. INTRODUCTION 

Laser control of molecular rotation, alignment and ori- 
entation has received significant attention in recent years 
(for a review, see e. g. 1, 2]). Interest in the field has in- 
creased, mainly due to the improved capabilities to ma- 
nipulate the characteristics of the laser pulses (such as 
time duration and temporal shape), which in turn leads 
to potential applications offered by controlling the angu- 
lar distribution of molecules. Since the typical rotational 
time scale is 'long' (~10 ps) compared to the typical short 
pulse duration (^50 fs), effective rotational control and 
manipulation are in reach. During the last decade, coher- 
ent rotational dynamics of pulse-excited molecules was 
studied % ^ , and multiple pulse sequences giving rise to 
the enhanced alignment were suggested [a-Qi ^^^ real- 
ized experimentally [SMllj. Further manipulations, such 
as optical molecular centrifuge and alignment-dependent 
strong field ionization of molecules, were demonstrated 
[12l |l3|. Selective rotational excitation in bimolecular 
mixtures was suggested and demonstrated in the mix- 
tures of molecular isotopes [1J| and molecular spin iso- 
mers [iSl Il6| . These new methods for manipulation of 
molecular rotation can also be used to modify the mo- 
tion of molecules in inhomogeneous fields, such as focused 
laser beams |17H19| . or static electric [20] and magnetic 
[lii fields. 

A molecule near a solid surface can be also considered 
as a particle in a complicated inhomogeneous field. Al- 
though the potential energy of a molecule near a solid 
surface is quite complicated in all its details, there are 
cases when simple potential models can be used. We 
treat below the molecule as a rigid rotor [23, while the 
surface is considered fiat and is described by a hard cube 
model [23, 2^. This model was used and gave quali- 
tatively correct results for N2 or NO molecules incident 
with thermal velocities on close-packed surfaces, such as 
Ag(lll) (see reference t25j] and references therein). 

Modification of the molecule-surface scattering and 
molecule-surface reactions by external fields of different 



nature is a long-standing research problem. In particu- 
lar, the effect of molecular orientation by a static electric 
hexapole field on the scattering process was investigated 
in detail |2^] . Laser control of the gas-surface scattering 
was achieved using multiphoton ionization of the imping- 
ing molecules by long laser pulses of variable polarization 
|2j|, and possibility of controlling molecular adsorption 
on solid surfaces using ultrashort laser pulses was dis- 
cussed J28|. 

In this paper, we investigate the prospects of modifying 
and controlling the process of molecular scattering from 
solid surfaces by using ultrashort laser pulses that align 
molecules before they hit the surface. 

The paper is organized as follows. The molecule- 
surface collision model is presented in Sec. II. Next, in 
Sec. Ill, the model for the interaction of a molecule with 
an ultrashort laser pulse is briefly explained. Then, we 
suggest two schemes for laser control of the molecular 
scattering process. In the first one, explained in Sec. IV, 
a single laser pulse is applied to a molecular beam in or- 
der to cause the molecules to rotate with a chosen sense 
of rotation (clockwise or counterclockwise) after the scat- 
tering from the surface. In the second suggested scheme, 
explained in Sec. V, two properly delayed laser pulses 
are applied to a molecular beam composed of two chem- 
ically close molecular species (isotopes, or nuclear spin 
isomers). As a result of the double pulse excitation, the 
subspecular scattering angles become enriched in one of 
these species after the scattering from the surface. In 
Sec. VI we summarize and conclude. 



II. MOLECULE-SURFACE SCATTERING 
MODEL 

In this paper, we use a model in which molecule is 
treated as a rigid dumbbell 24]. This dumbbell collides 
with a fiat frictionless hard cube, that represents one of 
the surface atoms [2J|. We assume that the cube has 
some velocity that is distributed according to the surface 



temperature. This hard cube model provides a simple 
way of adding surface phonons to the molecule-surface 
collision process. For the sake of simplicity, we also as- 
sume that the cube is much heavier than the molecule, 
so that its velocity does not change as a result of the 
collision. In this case, by moving to the frame attached 
to the cube, one reduces the problem to the molecular 
collision with a motionless hard wall. In the moving 
coordinate system, the molecular total energy (trans- 
lational-|-rotational) is conserved, but it can be redis- 
tributed between these two parts as a result of the colli- 
sion. 

What about the translational linear momentum of the 
molecule in the moving coordinate system? The compo- 
nent perpendicular to the surface is not conserved, be- 
cause the surface exerts forces on the molecule in this di- 
rection during the collision. On the other hand, there are 
no forces applied in the direction parallel to the friction- 
less surface, and therefore the linear momentum parallel 
to the surface is conserved. This is the reason why we can 
simplify the problem and consider a colliding molecule 
that has only a translational velocity component perpen- 
dicular to the surface. Notice that this cannot be done if 
the surface is corrugated. 

Using energy and angular momentum conservation 
laws (as explained below), we find analytic expressions 
for the translational and the rotational velocities of the 
dumbbell molecule after the collision. These velocities 
depend on the velocities before the collision, and on the 
angle between the dumbbell and the surface of the cube 
at the moment of collision. Finally, we transform the 
velocities back to the laboratory coordinate frame. 

In the next subsections, we treat a simple case of a 
homonuclear molecule rotating in a plane and colliding 
with a heavy hard cube. Next, we extend the treatment 
to a 3-dimensional rotation of a heteronuclear molecule. 
In the last subsection, we consider the effects of the sur- 
face cube vibration on the collision. 



The tw^o-dimensional collision of a homonuclear 
diatomic molecule 



We treat a homonuclear diatomic molecule as a mass- 
less stick of length re, with two atoms, each of mass m, 
attached to its ends. To describe the molecular motion, 
we define the Z-axis perpendicular to the surface, see 
Fig. [TJ The angle between the molecular axis and the 
Z-axis is 6, and it belongs to the range of [0,27r]. The 
translational velocity of the center of mass is denoted by 
V. The linear rotational velocity is denoted by v, and 
is equal to r'ew/2, where uj is the angular velocity of the 
molecular rotation. The velocities V and v are defined 
in the coordinate system moving with the velocity of the 
surface cube. 




FIG. 1: A homonuclear diatomic molecule hitting a flat hard 
surface. This is a two-dimensional model, where the molecule 
rotates only in the plane of the figure. The translational ve- 
locity of the center of mass is V, the rotational velocity of 
each one of the atoms, in the coordinate system moving with 
the center of mass, is v, and the angle between the surface 
normal (the Z-axis) and the molecular axis is 6. 



Energy conservation. The total energy conservation 
for the molecule is: 



i(2m)V^/ + ^lojf = i(2m)\f + llco 



f 



(1) 



where / = mrl/2 is the moment of inertia, and the sub- 
scripts i and / denote the velocities before and after the 
collision, respectively. Simplifying the above expression 
we obtain: 



V/ 



= v 



'f ' 



{V,-Vf){V, + Vf) = {vf-v,){vf+v.) 



(2) 



(3) 



Angular momentum conservation. The angular 
momentum of this system of two particles depends on 
the choice of the coordinate system. We choose it such 
that at the moment of collision the origin is at the po- 
sition of the colliding atom. Thus, the other atom is 
the only one contributing to the angular momentum of 
the system. Moreover, the torque exerted by the wall is 
zero for this choice of the coordinate system. Therefore, 
for this choice of the coordinate system, the angular mo- 
mentum is conserved during the collision. Adding the 
velocities V and v, as defined in Fig. [U and assuming 9 
being between and 7r/2, we equate the magnitude of 
the angular momentum before and after the collision: 

mfeivi — Vi sin 9) ~ mre{vf — Vf sin 9) , (4) 



{V,-Vf)sin9 = v,^Vf . (5) 

Dividing ([3]) by (O and using the result together with 
Eq. (O, we obtain the expressions for Vf and Vf: 



Vf = 



—Vi cos^ 9 — 2vi sin ( 



1 



sm 



«/ 



Vi cos^ 6 — 2Vi sin 6 
1 + sin^ e 



(6) 



These equations should give the same results ior 9 ^ 9 + 
TT, because of the symmetry of the molecule. Because the 
sine function changes sign under this transformation, the 
correct equations for 9 in the first or the third quadrant, 
are: 



Vf 



-Ficos^ 9 — 2vi\ sint 



1 



sm 



«/ 



Vicos^ 9 — 2Vi\ sint 
1 + sin^ 9 



(7) 



for < 6* < 7r/2 and tt < 6* < 37r/2 



For the angle 9 in the second or the fourth quadrant, 
the transformation — > tt — should be done. This is 
equivalent to changing Vi and Vf to —Vi and —Vf. Finally, 
we obtain: 



Vf 



-V^icos^ 9 + 2vi\ sin( 
1 + sin^ 9 



_ V, cos2 6' + 2F,|sin6'| 
"^ ~ 1 + sin^ 9 ' 

for tt/2 <9 <TT and 3tt/2 <9 <2tt 



(8) 



Notice that Vi is always negative, and Vf is posi- 
tive/negative for a clockwise/counterclockwise rotation, 
respectively. 

The angles 9=-k/2 or 37r/2 are unique and are excluded 
from the equations above. The reason is that the treat- 
ment above is incorrect for this angle of incidence, be- 
cause both atoms hit the surface at the same time, and 
the law of conservation of angular momentum, as applied 
above, cannot be applied in this specific case. However, 
these angles of incidence can be treated as limiting cases 
of a double collision, and treating them provides no par- 
ticular problem, as explained below. 

Time evolution. The above analytical expressions 
for translational and rotational velocities describe a sin- 
gle collision of the molecule with the surface. However, 
additional collisions may occur as well. This scenario is 
most evident for molecules hitting the surface at the an- 
gle close to 7r/2. In the case of the additional collision, 
the equations ([7]) and dH) should be applied again. 

Instead of finding complicated conditions for multiple 
collisions, we turn to a simple numerical simulation. The 
equations of motion for the molecule, starting at some Zq 
and 9o with velocities Vi and Vi, before the first collision, 
are simple: 



Z{t) = Zo + V^t , 



(9) 



and 



9{t) = 00 + Wot = 6*0 + —Vit, mod {2tt) (10) 



The collision occurs when the distance between the face 
of the cube and the center of mass of the molecule is equal 
to 0.5re| cos6'(i)|. At this time moment, the velocities are 
transformed according to ([7]) or (|S]), and the evolution 
continues according to ([9]) and (|10p , using the new values 
of Zo, Vi, 9o and Vi. We assume here that the collision 
lasts for a very short period of time, such that during the 
collision the values of Z and 9 are practically constant. 



B. Extension to three-dimensional collision of a 
heteronuclear diatomic molecule 

We describe the heteronuclear molecule as composed 
of two different atoms of masses mi and TO2, where 
TOi > 7712. The atoms are connected by a massless rod 
of length Te, as shown in Fig. [2][a). The center of mass 
of the molecule is closer to the heavy atom, and is at the 
distance oi re/{fi + l) from it, according to the definition 
of the center of mass. Here /x was defined as mi /?7i2. The 
distance between the center of mass and the light atom 
is, accordingly, re/x/(/L(,-|- 1). 

The molecule is characterized by the center of mass 
coordinate Z, and the direction of the molecular axis. 



described by a vector r 



'—re{x,y,z), where {x,y,z) 



is a unit vector. This vector points from the molecular 
center of mass to the light atom of mass TO2, see Fig. 
[2ja). The corresponding velocity is given by v = dv/dt. 
Therefore, the linear velocity (in the center of mass co- 
ordinate system) of the light atom is v, and that of the 
heavy atom is — v/^. 

After a derivation, similar to the one for the 2-D rota- 
tion presented above, we obtain expressions connecting 
velocities before the collision, V^ and v.;, to the veloci- 
ties after the collision, Vy and v/. The molecule, at the 
moment of collision, is oriented at r^, see Fig. [2{b). The 
details of the derivation are given in the Appendix at the 
end of the paper. 

For a free-rotating heteronuclear molecule, the end of 
the vector r traces a circle with a radius equal to — ^rfJ^e- 
The orientation of the circle in space is defined by the 
initial orientation ro and the initial velocity vq. The 
orientation of the molecule at time t is thus given by: 

r(t) = ro cos (wot) H — -rg— sin(uot) , (11) 

/i + 1 Vq 

where Vq ~ |vo|. Taking the derivative, the velocity is: 



v(t) == -t;orosin(uot) + 



A* 



Ai 



1 



reVo cos {vot) 



fl2) 



Finally, we emphasize that the molecular collision with 
the surface of the cube occurs when the molecule is ori- 
ented with Zi > Q and the distance between the face of 
the cube and the center of mass of the molecule is equal 
to re|cos0|/(/i -1-1), or when the molecule is oriented 
with Zi < Q and the distance between the face of the 
cube and the center of mass of the molecule is equal to 
re|cos6'|/x/(^-|- 1). 




Vf-Vi 




(a) 



(b) 



FIG. 2: (a) A heteronuclear diatomic molecule composed of 
atoms of masses mi and m2 (where mi > 7712) connected by 
a massless rod of length Vg. The translational velocity of the 
molecular center of mass is V, the molecular orientation in 
space is represented by the vector r, and the corresponding 
rotational velocity of the light atom is v. (b) A vector dia- 
gram corresponding to Eq. (|A.3|) . where molecular orientation 
vector at the moment of collision is denoted by r^. 



C. Vibration of surface atoms - simple inclusion of 
phonons 

The hard cube model [23, U^ [23 provides a simple 
way for including surface atom vibration into the colli- 
sion process. The collision of a molecule vifith a hard 
wall of infinite mass is replaced by a collision with a hard 
cube of a finite, but a large mass M moving with veloc- 
ity U. For simplicity, we assume that the cube oscillates 
in a hard box of a finite size, which is a free parameter 
of the model. In the model, the cube oscillates in the 
direction perpendicular to the surface plane. The reason 
is that for a fiat and frictionless cube only the vertical 
velocity component can transfer energy to the imping- 
ing molecule. The cube moves with a constant speed, 
while the velocity reverses its direction at the ends of the 
hard box. The one dimensional velocity of the hard cube 
is random and is distributed according to a Boltzmann 
distribution: 



f{U)^ 



M 



27rfc_BTsurf 



exp 






(13) 



where Tgurf is the temperature of the solid surface. When 
treating the collision process, we assume that the mass 
of the surface atom is much larger than the mass of the 
molecule, i. e. Af ^ ttii -|- TO2, so that the collision does 
not change significantly the velocity of the hard cube. 
Only the velocities of the molecule change, and its total 
energy in the laboratory coordinate frame may increase 
or decrease. 

The previous treatment of the collision in the coor- 
dinate frame moving with the cube is now easily incor- 
porated into the model. The collision condition (for a 
homonuclear molecule) can be expressed as 

Z{t)-ZMit)^^r,\cose{t)\, (14) 



where Zm (t) is the time-dependent position of the surface 
of the hard cube. According to the assumptions above, 
Z]\j{t) is a simple "zigzag" function (triangle wave) with 
an amplitude that is a free parameter, which we take 
equal to the molecular bond length r^, and a frequency 
determined by the cube speed \U\. The translational ve- 
locity of the molecule is transformed before the collision 
according to Vi ^ Vi — Ucoi, where Ucoi is the velocity 
of the hard cube at the moment of collision. This is a 
usual Galilean transformation to the coordinate frame 
moving with the hard cube. Similarly, after the collision 
the translational velocity of the molecule is transformed 
back by Vf ^ Vf + Ucoi- The rotational velocity remains 
unchanged under this Galilean transformation. 



III. INTERACTION OF THE MOLECULE 
WITH AN ULTRASHORT LASER PULSE 

Here we briefiy summarize the results of the classical 
model describing the interaction of the diatomic rigid 
molecule with a nonresonant ultrashort laser pulse, in the 
impulsive approximation. A more detailed description 
may be found in [30|. 

The potential energy of the laser pulse interacting with 
the induced molecular dipole is given by: 

V{9,ip,t)^~-j£^{t){Aacos^ l3 + aj_) , (15) 

where Aa — an — a± is the difference between the polar- 
izability along the molecular axis and the one perpen- 
dicular to it, £{t) is the envelope of the electric field 
of the linearly polarized laser pulse, and /? — j3{9, (p) is 
the angle between the molecular axis and the direction 
of polarization of the pulse. Here 0, ip are the polar and 
the azimuthal angles characterizing the orientation of the 
molecular axis, respectively. We assume that the pulse 
duration is very short compared to the rotational pe- 
riod, so that the pulse can be described in the impulsive 
((5-kick) approximation. We define the dimensionless in- 
teraction strength P, which characterizes the pulse, as 



„ Aa ; „T, s , 



(16) 



We consider the action of a pulse linearly polarized 
along some arbitrary unit vector p, and determine the 
vector of the resulting velocity change Av for a molecule 
oriented along some direction Tq. The norm |Av| can 
be found by integrating Newton's equations of motion 
for a pulse polarized along the z-axis. It is equal to 
j|Psin2/3o|, where /3o is the angle between the polariza- 
tion direction of the pulse p and the orientation direction 
of the molecule ro , and I is the moment of inertia of the 
molecule. Notice that Av is always perpendicular to Tq. 
Also, Av is directed parallel or antiparallel to the vector 
component of p perpendicular to ro, which is equal to 



p — -^^^ cos/3oro. As a result, we arrive at: 



Av ^ —^ cos/3o ( , , 7-eP - cos/jpro ) • (17) 



laser pulse 



^i + l 



IV. "MOLECULAR PROPELLER" INDUCED 

BY LASER ALIGNMENT AND COLLISION 

WITH THE SURFACE 

In this section we explore a way of inducing unidirec- 
tional molecular rotation by a single laser pulse and a 
single surface scattering event. The idea is inspired by a 
recent scheme that was proposed in [3l|, JS^] and realized 
experimentally in [33]. In these papers, two time-delayed 
and cross-polarized laser pulses were used in order to in- 
duce molecular rotation of a preferred sense. Here we 
achieve a similar goal by replacing the second laser pulse 
by the process of molecular scattering from a solid sur- 
face. 

Before presenting our new scheme, we summarize 
shortly the pure optical one described in |30l - l33{ . We 
start with a gas of diatomic molecules in free space. The 
first ultrashort laser pulse, linearly polarized along the z- 
axis, induces coherent molecular rotation that continues 
after the end of the pulse. The molecules rotate under 
field-free conditions until they reach an aligned state, in 
which the molecular axis with the highest polarizability 
is confined to a narrow cone around the polarization di- 
rection of the first pulse. The second short laser pulse 
is applied at the moment of the best alignment, and at 
angle with respect to the first pulse. As a result, the 
aligned molecular ensemble experiences a torque causing 
molecular rotation in the plane defined by the two polar- 
ization vectors. The rotational velocity delivered to a lin- 
ear molecule is maximal when the laser pulse is polarized 
at 45 degrees with respect to the molecular axis of the 
highest polarizability, as can be seen from Eq. ([T71) . This 
defines the optimal angle between the laser pulses. The 
direction of the excited rotation (clockwise or counter- 
clockwise) is determined by the sign of the relative angle 
(±45 degrees) between the first and the second pulse in 
the polarization plane. This double pulse scheme was 
termed "molecular propeller" , as it resembles the action 
needed to ignite a rotation of a plane propeller. 

In the current scheme, we start from a monoenergetic 
molecular beam of diatomic molecules flying towards a 
flat surface at the incidence angle of 45°. Before hitting 
the surface, the molecules are aligned by a laser pulse po- 
larized at -1-45° to the surface. When colliding with the 
surface, the aligned molecules receive a "kick" from it, 
and scatter with rotation in a specific direction. Chang- 
ing the polarization angle of the laser with respect to the 
surface to —45°, the sense of rotation of the scattered 
molecules can be inverted. 

In order to analyze properly the scattering of a molec- 
ular beam, we need to account for the spatial spread of 
the molecules inside the beam. For typical experimen- 




FIG. 3: (Color online) A molecular beam impinges on a hard- 
wall surface at angle of 45°. Before hitting the surface, the 
molecules are "kicked" by an ultrashort laser pulse polarized 
as shown by the double-headed (red) arrow. This pulse gen- 
erates the time-averaged angular distribution of the molecu- 
lar orientation in the form of the "cigar" as shown. These 
aligned molecules are preferentially rotating clockwise after 
hitting the surface. 



tal conditions, this spread is of the order of 1 mm in the 
direction perpendicular to the direction of the molecular 
beam. This spread is determined by the diameter of the 
collimating aperture in the experimental system [3J]. A 
similar spread can be observed in the direction parallel 
to the direction of the beam propagation, if one uses the 
pulsed molecular beam technique 35] . We ignore in our 
calculation the spread in the perpendicular direction and 
concentrate on the parallel one. 

After the molecules are "kicked" by the laser pulse, 
they start rotating in a concerted way, while continuing 
to approach the surface. The angular distribution of the 
"kicked" molecules, after averaging over a long time pe- 
riod, is elongated along the pulse polarization direction, 
as was shown in [30|, and as is depicted schematically 
in Fig. [3] by the "cigar" -shaped distribution. The result 
of this time averaging is the same as of the averaging 
along the distance parallel to the direction of molecu- 
lar propagation, because time is related to distance by 
Z — Vot/^/2 for all the molecules. From the above we 
conclude that on average the molecules approaching the 
surface are aligned along the polarization direction of the 
ultrashort laser pulse. We stress, that this is true classi- 
cally, as well as quantum mechanically. If the molecules 
are aligned at the angle of 45 degrees with respect to the 
surface, they have a preferred orientation while colliding 
with the surface. These molecules receive a "kick" from 
the surface, which leads to the preferred sense of rotation, 
for example, a clockwise rotation, as in Fig. |31 If we plot 
the distribution of the components of angular momentum 
of the scattered molecules, we expect to see an asymme- 
try that is correlated with the polarization direction of 
the exciting pulse. 

In the following, we apply the above scheme to a molec- 
ular beam of N2 molecules hitting a flat hard-cube sur- 
face with an angle of incidence of 45°. The molecules 
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FIG. 4: The distribution of the angular momentum compo- 
nents Jx , Jy and Jz of the molecules scattered from a hard- 
cube surface are plotted in the three panels. These molecules 
are manipulated with an ultrashort laser pulse before hitting 
the surface, as shown in Fig. [S] As expected, the distribution 
of Jy is asymmetric, and there are more molecules with a pos- 
itive Jy. The average value of Jy is (Jy) — 4.2ft. The calcula- 
tion is done with the hard-cube model for nitrogen molecules 
(atomic mass of 28a.u.) incident with translational velocity 
of 350 ^2- and rotational temperature of 1 K. The laser pulse 

sec ^ ^ 

strength is P = 10, and the polarization angle is 7 = 45°. 
The surface is composed of silver atoms at room temperature 
(M = 108 a.u. and T^u^t = 300 K). 



FIG. 5: (Color online) The percentage of the laser-excited 
molecules scattered to the subspecular angles as a function 
of the polarization angle of the pulse 7 (solid blue line). The 
dashed green line represents the average induced angular mo- 
mentum component (Jy) as a function of 7. The parameters 
are the same as in Fig. [l] 



of the pulse in the xz-plane. We denote this angle by 7, 
so that a pulse polarized along the z-axis corresponds to 
7 = 0. It is seen, as expected, that the largest induced 
(Jy) is obtained for a polarization angle close to 45°. 



move with the initial velocity of Vq = 350 -^ , and with 
the rotational temperature of 1 K, typical for molecular 
beam experiments. They receive a kick with a strength of 
P = 10 from an ultrashort laser pulse, before hitting the 
surface. We choose the surface to be Ag(lll), with the 
appropriate mass of the representative hard cube, and 
the surface temperature is taken to be 300 K. 

We choose the polarization direction of the exciting 
pulse to be p = (1,0, l)/\/2, that is, at 45° to the sur- 
face in the xz-plane, as shown in Fig.[3l Fig.|4]shows the 
distribution functions of the angular momentum compo- 
nents Jx, Jy and Jz after the scattering. We see that the 
distributions of J^ and Jz are symmetric around zero, as 
expected from the symmetry of the pulse. However, most 
of the molecules are scattered with positive Jy, and its 
average value is (Jy) = A.2h. This average is close to the 
value of of 5. 6/1, that may be estimated by using Eq. ([8]) 
for Vf, and the appropriate molecular constants m and 
Tg. This estimation considers a representative nonrotat- 
ing molecule oriented at 9 — 7r/4 and impinging on the 
surface with the incident velocity of Vi — 350 m/sec. 

Switching the polarization direction to p = 
(—1,0, l)/\/2 inverts the Jy distribution of the scattered 
molecules, and makes it peaking at negative values of Jy. 

In Fig. [5] we plot by the dashed (green) line the aver- 
age induced y-component of the angular momentum as a 
function of the angle of the linear polarization direction 



V. LASER-CONTROLLED SURFACE 

SCATTERING AND SEPARATION OF 

MOLECULAR MIXTURES 

In this section, we show that, in principle, one may use 
laser-controlled surface scattering for separating molec- 
ular beams consisting of several molecular species into 
individual components. The scheme seems to be appli- 
cable to different types of molecular species, such as iso- 
topes, or nuclear spin isomers. It takes advantage of the 
fact that rotationally excited and unexcited molecules 
have different scattering angle distribution after a colli- 
sion with the surface. 

As an example, we consider a molecular beam com- 
posed of a mixture of two nitrogen species. They can 
be two molecular isotopes, such as ^''N2 and ^^N2, or 
two nuclear spin isomers, such as ortho and para isomers 
of "'^^N2. It was shown in the past, both theoretically 
and experimentally, that two properly delayed ultrashort 
laser pulses may selectively align a preferred component 
of such a mixture, while lea ving the other one practically 
unexcited (see [ij and [l5l.[l6|). 

In order to explain this manipulation, a quantum de- 
scription of the kicked molecules should be used. After 
the molecules are kicked by a single laser pulse, they are 
transiently aligned, and shortly after that become ran- 
domly oriented again. However, because the quantum 




FIG. 6: (Color online) A mixture of two molecular species is 
manipulated by two properly delayed ultrashort laser pulses 
before hitting the surface. The species represented by dumb- 
bells with square ends (blue) are rotationally cold (their rota- 
tion is de-excited by the second laser pulse). The species rep- 
resented by dumbbells with circular ends (red) are rotation- 
ally hot (their rotation was further enhanced by the second 
laser pulse). After hitting the surface, the "blue" species are 
mostly scattered specularly, and towards angles larger than 
the specular angle. On the other hand, the "red" species 
mostly transfer their rotational energy to translational en- 
ergy, and are scattered to angles smaller than the specular 
angle. 



energy levels of the rotor are discrete, and because of 
the symmetry of a linear molecule, the molecular align- 
ment reappears later due to the phenomenon of quantum 
revival of the rotational wave packet. Generally, the dy- 
namics of the rotational wave- function repeats itself after 
a fixed time period, called the revival time, which is pro- 
portional to the moment of inertia of the molecule. 

It was shown in [ij] , that molecules aligned by the first 
laser pulse may become even more profoundly aligned if 
a second laser pulse is applied to them, which is delayed 
by an integer multiple of the revival time. However, if 
the time delay between the two pulses is close to a half- 
integer multiple of the revival time, the rotational energy 
given to the molecules by the first pulse is taken away by 
the second one, and the molecules practically return to 
the unexcited isotropic state they were in before the first 
pulse. 

Consider now a mixture of two isotopes mentioned 
above. Because of the mass difference of the isotopes, 
they have different revival times. After this mixture is 
kicked by the first pulse, it is possible to find delay times 
such that one of the isotopes evolved for an integer num- 
ber of the revival periods, while the second one completed 
a half-integer number of its own periods. If the second 
pulse is applied at one of these moments, then the first 
isotope will experience enhanced rotational alignment, 
while the second one will become isotropic and rotation- 
ally de-excited. 

A similar double pulse approach can be used for se- 
lective alignment in a mixture of nuclear spin isomers 
of "'^^N2, as explained in [1^]. In this case, both isomers 
have the same mass and, therefore, the same revival time. 
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FIG. 7: (Color online) The scattering angle distribution of 
the "square atom" (blue) species at the top and the "circular 
atom" (red) species at the bottom, according to the scheme at 
Fig. El The calculation is done using the hard-cube model for 
nitrogen molecules incident on the surface with translational 
velocity of 350m/sec and rotational temperature of 1 K. The 
distribution at the top panel is for the molecules that are not 
excited by the laser, in order to represent the species that are 
de-excited by the second laser pulse. The distribution at the 
bottom panel is for the molecules that are excited by a laser 
pulse with P = 10 and 7 = 0, to represent the species that are 
excited twice by a laser pulse with P = 5. In the upper panel 
of the figure, only 65% of the molecules are scattered to sub- 
specular angles, while at the lower panel, 84% are scattered 
to these angles. This means that almost 30% enrichment is 
achieved for the selectively excited isotope. The surface cube 
parameters correspond to silver atoms at room temperature 
(M = 108 a.u. and Tsurf = 300 K). 



However, kicking the molecules with a second pulse de- 
layed by a time close to a n -I- 1/4 (or a n -I- 3/4) multiple 
of the revival time (where n is a positive integer) provides 
further rotational excitation of one of the isomers, while 
de-exciting the other one. The reason for this effect is 
the entanglement between the molecular rotational and 
spin degrees of freedom, which is imposed by the Pauli 
principle (for details, see [l6j])- 

Assume now that such a selectively excited molecular 
beam hits a solid surface, as shown schematically in Fig. 
ini For a moment, we also assume that the surface is 
not vibrating, i.e. stays frozen at zero temperature. Ac- 
cording to Eq. ([T|), the rotationally hot (aligned) molecu- 
lar species, represented by dumbbells with circular ends 
(red) in Fig. [SI have a high probability of transferring 
their rotational energy to translational energy via the 
collision. As a result, the perpendicular component of 
molecular velocity after the collision is larger than the 
one before the scattering (the parallel velocity component 
remains the same). In other words, these molecules have 
a high probability to be scattered to the angles smaller 
than the specular angle (i.e., to the "subspecular" an- 
gles). The second species, that are rotationally cold (and 



isotropically oriented), represented by dumbbells with 
square ends (blue) in Fig. [HI have a high probability to 
transfer their translational energy to the rotational en- 
ergy, and they are mainly scattered to angles larger than 
the specular angle. As a result of the scattering pro- 
cess, the molecular mixture that arrives to the region of 
subspecular angles is highly enriched with the "circular 
atom" (red) species. At finite temperature of the sur- 
face, the effect is somehow reduced. The reason is that 
even a non-rotating molecule can acquire translational 
energy from a vibrating surface atom moving towards it, 
and scatter to subspecular angles. However, for a strong 
enough laser pulse, a sizable effect is expected even for 
surfaces at room temperature. 

To estimate the magnitude of the effect at typical ex- 
perimental conditions, we consider scattering of a beam 
consisting of a 1:1 mixture of two nitrogen isotopes dis- 
cussed above. The beam has initial translational velocity 
of 350^, rotational temperature of IK, and it propa- 
gates at 45° with respect to the surface. The surface 
consists of silver atoms, and the surface temperature is 
300 K. Before hitting the surface, the molecules are ex- 
cited by a pair of laser pulses of P = 5 that are polarized 
along p = (0, 0, 1) direction. The timing between the 
pulses is chosen such that one of the species remains ro- 
tationally unexcited, while the other one experiences an 
efficient kick of P = 10 from the double pulse. In Fig. 
[7] we plot the distribution of the scattering angle for the 
two molecular isotopes. The upper panel of the figure 
corresponds to the unexcited component of the molecu- 
lar beam, the lower panel corresponds to the selectively 
excited isotope. In both parts of the figure, the distri- 
butions are peaked around 25°. This peak comes from 
the translational energy delivered to the molecules by the 
thermally oscillating surface atoms. It is easy to under- 
stand the appearance of the maximum around 25° by 
using the following simple arguments. Consider a ball 
with velocity v, colliding with a heavy cube moving to- 
wards it with velocity U. This ball is reflected from the 
cube with velocity v + 2U. Treating the molecule as a 
ball, and keeping in mind that it hits the cube at angle of 

45°, we arrive at a scattering angle of arctan ( ""{J^^ j . 

Using a typical value of ^JkBTsmf/M = 150 ^ for the 
velocity C/, and the incident velocity v of 350 ^S- we esti- 

■^ ' ^ sec 

mate the scattering angle as 24° , which is consistent with 
the distribution maxima seen in Fig. [T] 

Analyzing the distribution functions of Fig. [71 we find 
that 84% of the rotationally hot species are scattered 
to the subspecular angles, while only 65% of the rota- 
tionally cold species are scattered to these angles, which 
means almost 30% enrichment in the selectively excited 
isotope. This figure is remarkable by itself, however it 
becomes even more impressive at lower surface tempera- 
ture. In Fig. |8] we plot the percentage of the subspecu- 
larly scattered molecules for the rotationally cold species 
(dashed blue line), and for the rotationally hot species 
(solid red line) as a function of the surface tempera- 
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FIG. 8: (Color online) The percentage of the subspecularly 
scattered molecules is plotted as a function of the surface tem- 
perature for the scattering of the rotationally hot molecules 
(solid red), and for the rotationally cold molecules (dashed 
blue) . The parameters chosen for the calculation here are the 
same as in Fig. [T] except the surface temperature. 



ture. We observe that for any temperature between zero 
and the room temperature, there are more laser-excited 
molecules scattered to the subspecular angles, than there 
are the unexcited ones. The effect is enhanced dramat- 
ically for surface temperature below 50 K, and the ra- 
tio between the two isotopes asymptotically tends to the 
impressive value of 70 at zero surface temperature! This 
suggests that cooling the solid surface increases the effect 
to a large extent. 

We also explored the dependence of scattering on the 
laser polarization direction 7 defined at the end of Section 
IV. In Fig. [5] we plot by the solid (blue) line the percent- 
age of the rotationally hot molecules scattered towards 
subspecular angles as a function of the angle 7. We find 
that the percentage varies slightly with the polarization 
direction and that the maximal percentage is obtained 
for pulses polarized at 7 = 0. 



VI. CONCLUSIONS 

We developed a simple classical model, based on the 
hard cube model, for the description of the molecule- 
surface scattering process. Using it and the classical 
model for interaction of molecules with ultrashort laser 
pulses, we suggested and investigated theoretically two 
possible schemes for modifying the process of molecule- 
surface scattering. 

In the first scheme, we proposed a way to exciting uni- 
directional rotation of molecules. First, the molecules in 
the molecular beam are aligned by a laser pulse in the di- 
rection at some angle with respect to the surface. Then, 
after the surface scattering, they were shown to rotate 



preferentially in a direction determined by the laser po- 
larization vector and the surface normal. This anisotropy 
of the angular momentum of the scattered molecules can 
be detected, for example, by using the REMPI spec- 
troscopy |33l.l36j. 

In the second scheme, we suggested exciting a molecu- 
lar beam consisting of two molecular species by two prop- 
erly delayed laser pulses, in order to provide selective ro- 
tational excitation of one of the species (similar to [ij] 
and fl6J)- We have shown, that these two species are 
scattered differently from a solid surface, in a way that 
allows for the enrichment of the scattered beam in one 
of these species for subspecular scattering angles. This 
result is potentially interesting for the analysis and sepa- 
ration of molecular mixtures of different kinds, including 
isotopes, and molecular nuclear spin isomers. 
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Appendix: Three-dimensional scattering of 
heteronuclear diatomic molecules from a hard wall 



Energy conservation. The conservation of the total 
(translational and rotational) energy before and after the 
collision is, similar to Eq. ([T|), 
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(A.l) 



Angular momentum conservation. The angular 
momentum conservation in the coordinate system cen- 
tered at the colliding atom should be expressed now in 
vector form. The molecule hits the surface when oriented 
at some r^ = -^re{xi,yi,Zi). The collision is very fast, 
such that the velocities change, but the orientation re- 
mains the same r/ = r^ in the course of collision. 

We start with the case when the molecule hits the sur- 
face with the heavier atom of mass 7771, or with Zi > 0. 
The angular momentum conservation gives, similar to 



Eq. Q): 



7772- 



M 



rj X (V, + Vj 



1712^^— -r,x (V/+V/) 



or: 



r^ X (V/-Vi+v/-Vi) = 



(A.2) 



(A.3) 



It follows from the last equation that the vector in the 
parentheses should be parallel (or antiparallel) to r^ . The 
vector Yf — V^ is perpendicular to the surface, because 
we only treat the velocity component perpendicular to 
the surface (the hard cube is flat and frictionless) . The 
vector v/ — v^ should be perpendicular to r^, because 
both Vi and v/ are perpendicular to r^. All this leads 
to the vector diagram in Fig. [2jb) . It follows from the 
diagram that: 



|v/ 



|Vf- 1^.1 sin ^ 



(A.4) 



an expression similar to Eq. ([5|) , where Vi and Vf denote 
the center-of-mass velocities with the appropriate sign. 



and sin 



yf (in the 3-D case, the angle 9 varies 



between and tt, so that sin0 is always positive). By 
squaring the last equation, we obtain 



'f 



2v, -v/ = (!//- y,)sin2 



(A.5) 



The final velocity Vf can be substituted from Eq. (jA.ip . 
The dot product can be found from the following expres- 
sion: 

(v/ - Vj) • v., = Vi • Vf - vf = |v/ ~Vi\vi Cy^ ■ e^j-^^ , 

where e„. and evf-vi are the unit vectors in the direction 
of Vi and v/ — v^, respectively. The unit vector ey^^y. 
can be easily found using Fig. [^fb): 



- (4 ■r,)ri 
1 - (e; • ri)^ 



(A.7) 



where e^ is the unit vector in the direction perpendicular 
to the surface (the direction of V/ — V^ in Fig. [2ljb)). 

Finally, combining Eqs. (jA.ip . (IA.4I) . (JA.SP and (|A.6p . 
we obtain for Vf, which can be positive or negative: 

_ (sin^ 9 ~ fi)Vi - 2\vi\ sin 6* {By. ■ Cy^-y^) 
^ ~ sin^ e + n ■ i ■ ) 



Using Eqs. (IA.4I) . (IA.7I) and (IA.8I) it is easy to find v/ as 
well. 

Now we analyze the case when the molecule hits the 
surface with the lighter atom of mass 7772, or with Zi < 
0. Here, the energy conservation is the same, but the 
angular momentum conservation is: 



-7771 TjX Vi V, 

fi \ ^i 



-7771 TiX Vf Vf 

M V M 

(A.9) 
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After a similar derivation, we find: 

(^sin^6'- l) Vi - 2|i;i|sin6' (e„, • Cyj-y.) 



Vf 



and 



/i sin^ 9 + 1 



|v/ - Vi| = ^|V/ - Vi\smt 



(A.IO) 
(A.ll) 



It can be easily checked that equations (jA.Sp . (IA.4I) . 
(JA.IOP and (|A.llj) reduce to equations dTj and dH) for 
/i = 1 and a two-dimensional rotation. 
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